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ABSTRACT 
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Context. We have obtained Herschel images at five wavelengths from 100 to 500 /im of a ~ 5.5 x 2.5 degree area centred 
on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time 
project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing 
on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured 
star formation. 

Aims. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared mor- 
phology, comparing it with features seen at other wavelengths. 

Methods. We use high-resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission fea- 
tures that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground 
. — Galactic cirrus. 

Results. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three 
arc-like structures extending out to ~ 21, ~ 26 and ~ 31 kpc respectively, in the south-western part of M31. The 
presence of these features, hosting ~ 2.2 x 10 6 M of dust, is safely confirmed by their detection in HI maps. Overall, 
we estimate a total dust mass of ~ 5.8 x 10 7 Mq, about 78% of which is contained in the two main ring-like structures 
at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a 
function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 
in the nucleus to ~ 275 in the outermost region. 

Conclusions. Dust in M31 extends significantly beyond its optical radius (~ 21 kpc) and what was previously mapped 
in the far-infrared. An annular-like segment, located approximately at R25, is clearly detected on both sides of the 
galaxy, and two other similar annular structures are undoubtedly detected on the south-west side even further out. 

Key words. Galaxies: individual: M31; Galaxies: ISM; Infrared: ISM 
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1. Introduction 

An overall view of the interstellar medium of the Milky 
Way is hampered by the huge column densities that any 
line of sight has to go through. Andromeda is the closest 
galaxy offering broadly similar characteristics and an un- 
encumbered view, thus probably being the best target to 
test our knowledge of the physical processes that govern the 
formation and evolution of massive spiral galaxies. In fact, 
high quality data of M31 ar e available at a wea lth of diffcr- 
ent spectral ranges: X -rays ( Pietsch et al. 2005 1 , ultraviolet 



nent would have a SE P that will peak at wavelengths l onger 



than 200 fim (Haas et al. 1998 Montalto et al. 2009). Due 



dThilker et al.[[2005|, optical (e.g. SDSS, Ha 



Devereux 



et al.||1994| or the PAndAS survey, see e. g. |McConnachie et 
al.||200 8[), near and far-infrared (FIR) Paas et al.[ |1998 ' 

Beaton et al. 



2007 
2006 



2009 



Krause et al. 
and 



Barmby et al.[ [20061 JGordon et al.j [2006 



radio (e.g. 
. Despite this, the observationa 



m prep.[), millimeter ([ Nieten et al. 
Thilker et al. 



2004 Braun et al. 



to the wavelength limitations of the data, no constraints 
could be put on its presence and, even less, its characteris- 
tics. The estimates of the dust mass were therefore believed 
to be only lower limits to the real value. With Herschel data 
we are in the position of setting more stringent constraints 
both to the presence and spatial distribution of a very cold 
dust component, and its emissivity. 

The data we exploit in our survey allow us to also study 
the dust morphology in the very outer regions of the galaxy. 
The far-infrared morphology of Andromeda has been the 
topic of several studies exploiting, over the years, the avail- 
able infrared facilities such as the Infra Red Astronomical 



Satellite (IRAS, Habing et al. 1984 Xu fc Helou 1996), 



the Infrared Spatial Observatory (IS07Haas et~ 



coverage in the sub 
millimeter domain was until now limited by the absence 
of data longward of 170 /im, except for the extremely low 
resolution (with a beam size larger than 40') data from the 
Diff use Infrared Backgroun d Experiment (DIRBE) at 240 



and Spitzer ( Barmby et al. 2006 Gordon et al.l 2006 ) each 



1998) 



/im ( |Odenwald et al.[ [l998) 



facility improving in both sensitivity and spatial resolution. 
Evidence for a complete dusty ring at a radius of ~ 10 kpc 
came from the very first IRAS observations, while follow- 
ing studies exploiting higher resolution were able to detect 
details such as a discontinuity, most likely cause d by an in 



This situation has left M31 unexplored at far-infrared 
and sub-mm wavelengths in adequate spatial resolution. 
The adve nt of the Herschel Space Observatory ( Pilbratt et 
al. 2010), with its two instruments PACS (Photodetector 
Array Camera and Spectrometer) (Poglitsch et al. 2010) 



ter action with the satellite galaxy M32 (see, e.g. Gordon et 



al. 2006). A fa inter ring-like stru cture at ~ 15 kp c was aT 
ready noted by|Haas et al. ( 1998 1 and confirmed by Gordon 



and SPIRE (Spectr al and Photometric Imaging Receiver) 



(Griffin et al. 2010), allows us to fill this observational gap 



et al. ( 2006 ), but nothing more could be said on larger radii 
because of the size limitations of the maps. 

While dust located beyond the regions sampled by such 
observations was indeed detected indirectly by studying the 
characteristics of the stellar populations in the extreme out- 



and complete tETview of our neighbour galaxy down to skirtg of M31 (e g | C uillandre et al 



sub-mm wavelengths. The other important advantage of 
Herschel is the gain in terms of spatial resolution which, 
at 160 /im, is about 2.5 times higher than that of the 
Multiband Imaging Photometer for Spitzer (MIPS) at simi- 
lar wavelengths, making it possible to resolve physical scales 
of about 50 pc at the distance of Andromeda. 

On the one hand, M31 is considered in many respects 
to be very similar to the MW, sharing some characteristics 
and similarities such as the Hubble type, the luminosity, 
and gas content in the respective disks. The fact that it is 
so close makes it not only the ideal candidate for a compar- 
ison with similar studies and ongoing surveys tak ing place 
with Hersch el on our own Galaxy (such as Hi-Gal, Molinari] 
et al. 2010 ), but also the ideal place for studying the inter- 
stellar medium of a spiral galaxy. Hi nting at some com mon 



2001 



2012), it had never been possible to study its properties 



Bernard et al. 



aspects between the MW and M31, Yin et al. d2009b find 



very similar gas mass fraction profiles (similar scale lengths 
but different normalizations) and abundance gradients in 
the two disks. 

On the other hand, some differences are evident between 
the two galaxies: not only does M31 have about twice the 
baryonic mass of the MW, it has also a larger disk, by a 
factor of about 2.4. The mass of gas contained in the re- 
spective disks, however, is remarkably similar. An overview 
of the main prop erties of the two galaxies can be found in 
Yin et al.| ( [2009| (see their Table 1). 

So far, all the studies of M31's far-infrared Spectral 
Energy Distribution (SED), mainly relying on data short- 
wards of ~ 160 /im, have speculated about the possible 
presence of a very cold component (~ 10 K). Such a compo- 



* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with important participation from NASA. 



and quantify its presence before. 

In this work we describe Herschel PACS and SPIRE 
observations of M31, at 100, 160, 250, 350 and 500 /mi, ex- 
tending out to about 34 kpc (~ 2.5°) from its centre, pre- 
senting the most complete far-infrared/submillimeter sur- 
vey of Andromeda. We describe the observations and the 
data reduction steps in § [2j while in § [3] we give a gen- 
eral overview of the ancillary dataset that we exploit in our 
project, together with a comparison to Herschel maps. In § 
[4]we use high and low-resolution Hi maps, entirely covering 
our fields, to disentangle the dust dynamically connected to 
Andromeda from the foreground Galactic cirrus, and in §[5] 
we use these foreground-corrected maps to study the dust 
properties. Conclusions and remarks are given in § 6] 

Throughout this paper and in all works within the 
HELGA pr oject, we will assume a dist ance to Andromeda 
of 785 kpc ( |McConnachie et al.[ |2005|) , which turns into a 
linear scale of 3.8 pc/". As for the inclination and position 
angles, we are using i — 77° and PA= 38°, respectively. 

2. Observations and data reduction 

In this section we describe the observing strategy used to 
obtain the maps presented in this paper, and the data 
reduction techniques. Furthermore, given the greater dif- 
ficulty of removing instrumental artefacts in PACS data, 
some quality checks are performed for the latter. 



2.1. Observations strategy 

We observed a field of ~ 5.5° x 2.5° centered on M31, 
with both PACS and SPIRE in parallel mode, at a scan- 
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Fig. 1. Images in the five Herschel bands (from top to bottom: PACS 100 and 160 /xm, SPIRE 250, 3 50 a nd 500 fim). 
These are the original images, where no correction for the foreground cirrus was applied (see Sec. 4.1l. The green 
rectangles in the 250 /im map represent tentative detections of dust in the outskirts (see discussion in Sec. 4.2). The 
directional axes length corresponds to 15' (~ 3.42 kpc). o 
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ning speed of 60"/s. PACS data were taken at 100 and 
160 /im. Images at the five wavelengths are presented in 
Fig. [l] To obtain such a large map, we split our obser- 
vations in 2 parts, A and B, each observed with a cross- 
scan, so that the final map is a composition of four ob- 
servations. Each one covers a field about 3° x 2.5° wide, 
displaced approximatively along the major axis of the 
galaxy, and they overlap at the nucleus. The fields are cen- 
tered at RA=0h37m38.870s, DEC=+40dOQm51.30s {M 
and RA=0h48m01.670s DEC=+42d30m33.20s (eQ), and 
they both observe a stripe, parallel to the galaxy's minor 
axis, which is about 0.2° and ~ 0.3° wide in PACS and 
SPIRE maps, respectively. These observations were sched- 
uled on the 18th, 19th, 20th and 21st of December, 2010. 
Due to some issues with the star trackers, the last orthogo- 
nal scan of the northern portion of the map, was lost. The 
field was hence re-observed on the 23rd of January, 2011, so 
that it displays a slight rotation with respect to the nominal 
scan. 



2.2. PACS data reduction and quality check 

The basic PACS data reduction (up to level- 1 where pixel 
flagging, flux density conversion, and sky coordinate asso- 
ciation for each pixel of the detector arc applied) was per- 
formed in the Herschel Interactive Processing Environment 
(HIPE) v8, which includes the most recent responsivity fac- 
tors. At this stage, the PACS timelines are still dominated 
by brightness drifts caused by low-frequency noise. In or- 
der to subtract these drifts, remove glitches, and project 
the timelines on the final map, w e applied the IDL al go- 
rithm Scanamorphos (version 15, Roussel et al. 2011 ) to 
the level- 1 PACS timelines. To this end, it exploits the fact 
that each sky pixel is observed multiple times by several dif- 
ferent detectors, avoiding assumptions or modeling of the 
noise properties. The final PACS maps have a pixel scale of 
2" and 3" in the 100 (green) and 160 /im (red) band maps, 
respectively, and the FWHM of point sources is about 12.5" 
and 13.3" in the green and red channel, respectively, (see 
the "PACS photometer point spread function" report, by 
D. Lut^J. Note that, for high scan speed observations, the 
point spread function (PSF) is elongated in the scanning di- 
rection. This yields a spatial resolution on a physical scale 
of about 48.1 pc at 100 /zm and 50.7 pc at 160 /im.. 

In order to examine the accuracy of the PACS calibra- 
tion and of the flux output from Scanamorphos, we per- 
formed a pixel-to-pixel comparison between the PACS 160 
/xm and the Spitzer MIPS 160 /im images on the one hand, 



and b etween the PACS 100 /im and DIRBE ( Hauser et al. 
1998 1 at the same wavelength. We first convolved the PACs 



160 /im map (resp. PACS 100 /im map) to the MIPS 160 
/xm (resp. DIRBE 100 / im) resolution using the convolu- 



tion kernel provided by Aniano et al. (20111 (a gaussian 
kernel was used to smooth 100 /im data to the DIRBE 
resolution), and then resampled the maps on the same as- 
trometric grid using the hastrom task in IDL. Finally, we 
compared the flux in each pixel of the map. P ACS 100 /im 
and DIRBE 100 agree within 10%, as found in Meixner et 



1 OBSIDs 1342211294 and 1342211309 



2 OBSIDs 1342211319 and 1342213207. Data arc available at: 
http: //herschel . esac . esa. int/Science_Archive . shtml 



A November 3, 2010," version 1.01: herschel . esac . esa. int/ 
twiki/pub/Public/PacsCalibrationWeb/bolopsf vl . 01 .pdf 



al. (2012). As for the 160 /im band, we confirm the non- 
Tmear be havior of MIPS above 50 MJy/sr, which was also 
found by Stansberry et al. (20071 and is discussed in de- 
tail in [ Meixner et al | q2012} rBelow this threshold, PACS 
and MIPS agree within 5%7whereas the ratio PACS/MIPS 
is non-linear above that, leading to differences in bright- 
ness up to 30%, with MIPS underestimating the fluxes. 
Overall, the PACS extended calibration can be considered 
to be accurate to within 10%. A comparison of the total 
flux of M31 measured from PACS and MIPS maps at 160 
/im, indeed yields an agreement to within about 5%. We 
have also investigated diffe rent baseline subtrac tion meth- 
ods (HERITAGE-like, see |Meixner et al.[|2012| ) and found 
similar results. 

Hence, in the rest of this work, we will adopt a 10% 
value as uncertainty in PACS fluxes at both wavelengths. 

2.3. SPIRE data reduction 

The SPIRE data were processed up to Level-1 (i.e. to 
the level where the pointed photometer time-lines have 
been derived) with a custom driven pipeline script adapted 
from the official pipeline (P 'OF '5 -pipeline. py, dated 8 Jun 
2010) as provided by the SPIRE Instrument Control Centre 
(ICCjJ 

This Jython script was run in HIPE coming with the 
continuous integration build number 4.0.1367. In terms of 
the SPIRE scan map pipeline up to Level-1, this was very 
similar to the Herschel Common Science System/Standard 
Product Generation v5 since we already had patched the 
new parts of v5 into our old v4. So we used the new 
flux calibration product of v5 where the SPIRE calibra- 
tion is finally based on Neptune data, as well as the 
concurrentGlitchDeglitcher: the latter removes the ef- 
fect of often low level, but frequent glitches appearing in all 
bolometers of an array at the same time. 

The differences from the standard pipeline are that we 
used the sigmaKappaDeglitcher task instead of the ICC- 
default waveletDeglitcher and we did not run the de- 
fault temperatureDrif tCorrection and the residual, me- 
dian baseline sub traction. Instead we u se a custom method 
called BriGAdE ( |Smith et al.[ in prep. ), to remove the tem- 
perature drift arid bring all bolometers to the same level 
(equivalent to baseline removal). BriGAdE removes the 
temperature drift by fitting the thermistor timeline to the 
bolometer signal timeline assuming a linear relationship to 
fin d the baseline to subt ract (a similar process is presented 

In addition to BriGAdE's stan- 



Pascale et al. 



2010) 



dard source masking we also manually mask Andromeda 
from the timeline fitting and adjust the length of the lin- 
ear fit to every second scan-leg (including turn-around) to 
avoid issues interpolating the thermistor fitting process for 
the two observations where the satellite turns-around while 
pointing at the centre of Andromeda. For the two north- 
ern observations a gradient in the cirrus was found (as seen 
in previous IRAS observations) in the north-south direc- 
tion, however due to baseline subtraction this is not seen 
in the east-west observation. To compensate for this we use 
the IRAS map to identify a roughly constant level of cir- 



4 See "The SPIRE Analogue Signa l Chain and Photom eter 
Detector Data Proc essing Pipeline" 



Griffin et al. (20091 or 



Dowell et al. (20101 for a more detailed description of the 



pipeline and a list of the individual modules. 
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Fig. 2. U pper row, from left: the SPIRE 250 /im ima ge (this work), and far and near-ultraviolet images taken with 



GALEX ( |Thilker et alj |2005[ [Gil de Paz e t a l. [ |2007| ). Lower row, from left, the IRAS 100 fxm (MJy/sr) (shown as a 
com parison), the CO m ap (by Nieten et ai. 2006 in K km/s), and the Hi column density map (expressed in cm -2 ), 



from Braun et al. ( 2009 ) . All the images are set to the same spatial scale 



rus across the two observations and use these regions to 
perform our thermistor fitting. The northern and south- 
ern observations are then normalised to have the same zero 
level and mosaicked to create the final map. We have found 
this method improves the baseline subtraction significantly, 
especially in cases where there are strong temperature vari- 
ations during the observation. 

Our final maps were created using the naive mapper 
provided in the standard pipeline using pixel sizes of 6", 
8", and 12" at 250, 350, and 500 ^m, respectively. The 
non-standard pixel sizes were chosen in order for the pix- 
els to have approximatively the size of 1/3 of the beam. 
The FWHM of the SPIRE beams vary as a function of the 
pixel size and are 18.2", 24.5", and 36.0" at 250, 350, and 



500_um, for the pixel scale we have used (Swinyard et al. 
2010 1. In addition the 350 /im measurements are multipliec 



by 1.0067 to update our flux densities to the latest v7 cali- 
bration product. The calibration uncertainties are taken to 
be 5% correlated error between bands with a 2% random 
uncertainty, however it is recommended that these values 
are added linearly instead of in quadrature (see also SPIRE 
user manuaQ. 



5 http:// herschel.esac.esa.int/Docs/SPIRE/pdf/ 
spire_om.pdf 



3. Ancillary data 

In this section we first give a brief overview of the multi- 
wavelength dataset that is already available in the literature 
and that will be exploited by our project in future papers. 
We then present a comparison of Herschel data with other 
wide-field observations at various wavelengths, focusing in 
particular on the Hi datasets. 

In Fig. [2lwe show maps at UV, FIR (IRAS at 100 /im), 
millimeter " 2 CO(1-0)) and radio (Hi) wavelengths, com- 
pared to_SPl_RE^0 um (upp^r z le|tjar^._GALEX im- 



ages ( |Thilker et al.| 5 
PSF ranging from 4 , 
UV (FWHM), span r 



2005 
for 



Gil de Paz et al. 
'^ar-UV band, to 5' 



2007], with a 
for the Near- 



3° along the major axis, and they 
also include the satellite galaxy NGC205. In both maps the 
10 kpc (~ 0.75°) ring is clearly traced by UV emission, 
which extends out to the 15 kpc (~ 1.1°) ringlike struc- 
ture. The major difference from a morphological point of 
view, is found in correspondence of the "hole" w hich splits 
the ring {Gordon et all 120061 iBlock et al.[|2006|), which in 



UV images is peaking towards the outskirts. Otherwise UV 
and IR emission are well mixed across the galaxy. 

M31 was observed by Spitzer with bot h the InfraRcd 
Array Camera (IRAC) dBarmby et al||2006[ wit h a FWHM 
of 1.9"), and with MIPS {Gordon et al.[ |2006"1 PSFs from 
6" to 40"). The area imaged by the Spitzer survey was a 
^3x1° region, about the same size as that of the UV data. 
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CO data (Nieten et al. 20061, used to trace the molec- 
ular gas phase, have a spatial resolution of 23" and only 
cover a ~ 1.8° x 0.5° wide region, sampling the whole 10 
kpc ring. As expected, CO emission shows a remarkable 
morphological similarity with IR structures: faint nuclear 
emission with most of the radiation being emitted by the 
10 kpc ring. 

Finally, Hi maps, crucial to distinguish whether the dust 
belongs to M31 or to t he MW, are d i scussed in detail in 
Sectio n |3Tj both sets j| I liilkcr ct "iu~| |2004| |Braun et ~ 
2009) fully include the whole HELGA held (see e.g. Fig 



3.1. The Hi dataset 

Exploiting the velocity information carried by Hi observa- 
tions is a commonly used means to disentangle FIR emis- 
sion associated with Galactic cirrus, from that coming from 
extragalactic sources. Given the relatively low galactic lat- 
itude of Andromeda (b ~ —22°) contamination by cir- 
rus emission may be important in the field and has to 
be identified. Two datasets are used for this purpo se: a 
high-resolution cube, presented in Braun et al. fl2009|), and 
a low- resolution cube, presented in Thilker et al. ~ |2004 1 . 
Whenever possible, we used the high-resolution one. 

The high-resolution cube is derived in |Braun et al. 



( 2009 ) from mosaic Westerbork Synthesis Radio Telescope 
(WSRT) observations. The maximum angular resolution 
is 18" x 15", which corresponds to a spatial resolution of 
68 x 57 pc. The velocity coverage of the created maps 



ranges from -621 km s to -20 km s , and Braun et al. 



( 2009 ) report Galactic emission from -130 km s _i to +45 
km s _1 . However, we also make use of the low-resolution 



cube obtained by Thilker et al. ( 2004 ) from the Green Bank 
Telescope (GBT) observations, because (1) their velocity 
range (from -828 km s _1 to +227 km s _1 ) completely cov- 
ers the range of the foreground Galactic emission, and (2) 
foreground Galactic emission is expected to be angularly 
very extended, and GBT data are more sensitive to ex- 
tended emission than WSRT data. The GBT beam (9.1') 
corresponds to ~ 2.1 kpc at M31 distance. 



4. Results 

4.1. Foreground emission 

An attempt to detect diffuse dust structures in the outskirt 
of nea rby galaxies has been previously made by |Davies et al.| 
(2010), who carried out a detailed analysis of the infrared 



emission in the M81 group. They looked for an association 
of the filamentary emission, detected at SPIRE wavelengths 
in the proximity of M81, with similar structures detected in 
the optical. To identify the origin of this diffuse emission, 
they exploited the association of dust with atomic hydrogen 
and the kinematic information contained in the Hi dataset. 
They found that most of the diffuse far-infrared and optical 
emission located outside the clear extent of the disk of M81 
could be identified with Galactic cirrus. 

We have followed a similar approach: in order to dis- 
entangle pristine M31 emission from foreground Galactic 
(cirrus) emission, we have searched for identifications of 
the various IR structures visible in Fig. [5j with similar 
structures visible in the high-resolution Hi dataset. The 
latter has the advantage of carrying kinematic information: 



in principle, IR structures identified with Hi emission at 
non-Galactic velocities can be safely attributed to M31. 
However, the situation for M31 is somewhat unfortunate: 
its systemic velocity is approximately between -310 and - 
300 km s _1 and its maximal rotation velocity is of the order 
of 250-300 km s _1 (the inner rotation velocity, where the 
pea k rotation velocity is si tuated, is quite unc ertain, see 
e.g. IChemin et al.||2009| and |Corbelri et al ||2010[ ). This im- 
plies that for the receding (north-east) side of M31, where 
the velocity is close to km s _1 , disentangling Galactic and 
M31 emission is potentially problematic. 

To try to overcome this issue, we have derived maps 
of both M31 and the foreground Hi emission from the low 
resolution GBT observations. 

The foreground emission map, shown in the left panel of 
Fig. [3j was made in the appropriate range of velocities, from 
— 136 to +42 km s _1 , that includes the Galactic emission, 
under the assumption of negligible self-opacity over this in- 
terval. This foreground emission map has had regions of 
M31 emission replaced within each relevant velocity chan- 
nel, with a polygon containing a constant value related to 
the median value of the surrounding region at that same 
velocity. Th e polygons were defined, as described in Braun 



et al. ( 2009 1 , to contain all of the contiguously connected Hi 



emission m each velocity channel that could be attributed 
with confidence to M31. 

The foreground Hi image can then be used to determine 
appropriate scale factors to allow a correction of each of the 
FIR images. To achieve this, we have explored two differ- 
ent approaches. Firstly, using the Hi cirrus map, we have 
derived maps of the Galactic cirrus emission at Herschel 
wavelengths, by assuming a Galactic gas-to-dust ratio, a 
given dust temperature (we have tried a set of temperatures 
ranging from 10 to 15 K), and that the dust is emitting as 
a modified bla ck body with an emissivity value taken from 
Draine| (|2003|). Subtracting these images from SPIRE data 



to obtain cirrus-free maps, gave quite poor results, yielding 
heavy over-subtraction, with large portions of the maps dis- 
playing negative fluxes after this correction, mostly in the 
north-east region, the most confused area of the field. 

Given the potentially very large parameters space which 
is involved in this procedure (/3, dust temperature, gas-to- 
dust ratio, dust emissivity), making the solution to this 
problem highly degenerate, we have decided to try a differ- 
ent, more empirical, approach. We have looked for a pixel- 
by-pixel correlation between the FIR emission at SPIRE 
wavelengths and the Hi cirrus emission (Fig. [3l top panels) 
outside of M31 which we have derived as explained above, 
in maps that were convolved to a common angular reso- 
lution and matched in pixel size. Fig. [3] shows the scatter 
plots of the FIR emission versus Hi column density, to- 
gether with a linear fit to the points, at 250, 350 and 500 
/im, respectively. 

In these plots we note a group of data points displaying 
a Hi-FIR correlation with a clearly different slope with 
respect to our best fit, above Hi column densities of about 
10 21 cm -2 . Most of these points correspond to regions that 
are clearly associated with MW cirrus emission in the N-E 
side of Fig. [I] that is still be visible in the maps even after 
the correction (e.g. see Fig.|3j right side of the bottom row). 

The linear fit was performed in a limited range of Nhi 
and flux values, that was relevant for the foreground estima- 
tion. We have experimented extensively with other types of 
fitting, which included more of the FIR-bright regions, but 
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Fig. 3. Top row: on the left-hand panel the Hi galactic emission, as derived from GBT observations, is shown (see text 
for details). On the right-hand panel, M31 associated Hi emission, obtained as the difference between the GBT map 
integrated over all the velocities, and the pure MW cirrus emission. Bottom row: on the left-hand panel the original 
SPIRE 250 (im image, convolved to the GBT resolution and matched to the same pixelsize, compared to the cirrus 
subtracted image, at the same resolution (right-hand side). 



these always resulted in over-subtraction of the foreground 
in the most Northern portion of the FIR images. Cirrus 
maps at the SPIRE wavelengths are hence derived by mul- 
tiplying each pixel of the Hi foreground map by a value 
derived from the aforementioned correlations. A residual 



emission in all the FIR maps is still clearly visible in this 
part of the field (see e.g. Fig. [3] and Fig. [5|, having proba- 
bly very different dust /gas properties relative to the average 
Galactic foreground toward the M31 disk. Furthermore, we 
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Fig. 4. Pixel-by-pixel correlation between Hi column density, derived from the full Hi integral M31 image at low reso- 
lution, and the flux at 250, 350 and 500 /im from left to right, respectively. The red lines represent a linear fit to the 
data. 



stress that given the low spatial resolution of the Hi map, 
any cirrus structure smaller than ~ 9' will not be removed. 

The bottom panels of Fig. [3] show the comparison be- 
tween the original 250 /im map (left panel) and the cirrus- 
subtracted one (right panel; both images are at the GBT 
resolution) that has been obtained by applying the pixcl- 
by-pixel correlations mentioned above. The south-east edge 
of the latter image displays now some additional emission 
which was not present before the subtraction. This is an 
artifact, caused by the foreground subtraction, which only 
affects this particular region. In fact, this region has the 
lowest Hi column density (see upper-left panel of Fig. [3]) 
which, in the S^o^m vs Hi relation (Fig. [4|, lies in the 
negative side of the ordinate axis (this is due to slight im- 
perfections in the baseline subtraction). Hence, subtracting 
these negative values gives a positive contribution on this 
side of the map, but it is not affecting the emission from 
M31. 

The lines that produced the best overall corrections are 
overplotted in Fig. |4j The slopes of the lines for the 250, 
350 and 500 /im bands (in data where the beams have 
been matched to the Hi beam) are 4.50, 2.15 and 1.10 
Jy/beam/10 20 cm -2 , respectively. The relative FIR bright- 
nesses of the foreground in these bands are 1 : 0.48 : 0.24. 
For comparison, a modified blackbody with a dust tempera- 
ture of 15 K and and j3 of 2 would give relative brightnesses 
of 1 : 0.56 : 0.21, which is reasonably close to the mea- 
sured values and provides a physical justification for our 
method. We have used these correlations to derive, from the 
GBT Hi foreground image, FIR maps of the Galactic cir- 
rus emission at all SPIRE wavelengths, that we have then 
subtracted from the respective Herschel maps to correct 
for the presence of foreground dust. This correction only 
marginally affects M31 emission, yielding a pixel-to-pixel 
median difference of about 3% at all SPIRE bands, within a 
galacto-centric distance of ~ 20 kpc. A similar investigation 
on PACS data yielded no positive results, probably due to 
the lower sensitivity at these wavelengths. PACS maps are, 
hence, not corrected for the cirrus emission. 

In the rest of this work we use the foreground-corrected 
maps, unless stated otherwise. 
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Fig. 5. The Herschel 250 /im image processed for subtract- 
ing the foreground cirrus emission (see text), is shown. Dark 
green rectangles represent identifications of structures of 
which we investigate the association to M31. Blue contours 
represent the total column density for discr ete and diffuse 
M31 h igh-velocity Hi clouds, presented in Thilkcr et al. 
(2004) (Hi emission from the disk of Andromeda was re- 
moved), where we also look for the presence of dust asso- 
ciated to M31. Contours are 0.5, 1, 2, 5, 10, and 20 x 10 18 
atom cm -2 . The magenta star identifies Davies' Cloud. 



4.2. Far-Infrared structures 

Figure [I] presents images in the five Herschel bands. The 
morphology is characterized by a set of rings and arc-like 
structures whose prominence varies as a function of the 



1984 


Haas et al. 


1998 


Gordon et al. 



prominent structure at all wavelengths, while the inner- 
most regions become progressively less luminous at increas- 
ing wavelengths. Hints of a ring-shaped structure extend- 
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ing further out (~ 15 kpc) already seen both in ISO and 
Spitzer maps, can be noted at 100 /im, and its evidence be- 
comes clear in maps at all other wavelengths. Other even 
more external structures can only be seen in SPIRE images. 
Two other ring-like structures extending to the south-west 
are well detected at these frequencies, and they will be dis- 
cussed below. 

Images at all Herschel wavelengths are dominated by 
the 10 kpc ring which is clearly visible up to 500 /im. 
Furthermore, extending out to ~ 15 kpc from the centre, 
both PACS and SPIRE maps reveal a ring-shaped struc- 
ture, whose pre sence, at infrared wav elengths, had already 



been noticed by|Haas et al. (19981 and|Gordon et al. (2006) 



(see their Figs. 1). In the analysis performed within the iat- 
ter work, which exploits Spitzer images, it was already ques- 
tioned whether the far-infrared emission from Andromeda 
could extend to larger radii than those sampled by their 3° 
map. 

In Fig. [5] (but see also Fig. [T| possibly positive gas-dust 
associations with the disk and extended components are 
shown as rectangles. The north-east edge of M31 is con- 
taminated the most by foreground Galactic emission, and 
hence difficult to analyze, even after correcting for the fore- 
ground dust component. The two structures that we have 
indicated as C and D in Fig. [5] are spatially coincident 
with a structure that is probably associated with M31 and 
a structure that is probably foreground Galactic emission, 
respectively. This is illustrated in Fig. [6] For structure C, 
this is visible in channels —135.41 to —102.43 km s _1 (Hi 
emission associated with M31) and —52.96 to —44.72 km 
s _1 (Galactic foreground). For structure D, M31 emission 
is seen in channels —127.16 to —77.70 km s -1 , whereas in 
channels -69.45 to -44.72 km s" 1 the Hi emission from 
structure D is connected (spatially and in velocity) to both 
M31 emission and large-scale Galactic emission. 

The two structures that we have labelled as A and B 
are probably mostly due to Galactic foreground. Indeed, 
as Fig. [6] shows, there is a tiny amount of M31 Hi emis- 
sion at the position of structure A in channels —135.41 and 
—94.18 to —69.45, but there is relatively bright Galactic Hi 
emission in channels —52.96 to —19.98. Similarly, Hi at the 
position of structure B is seen in channels —52.96 to -44.72 
and —28.23 to —19.98. Because of the amount of extended 
Galactic Hi emission all over the map, it is probably fore- 
ground emission but, because of its position, it cannot be 
excluded that M31 contributes to some of the observed IR 
emission from that region as well. 

The unprecedented extent of our maps, together with 
the sensitivity of SPIRE, makes it possible to detect 3 fur- 
ther structures, extending to the south-west, that are most 
probably physically connected to M31. They appear to em- 
anate from the outermost ring, going out to a distance of 
~ 21, ~ 26 and ~ 31 kpc from the nucleus (see Fig. fl] [5] and 
[8j where they are named E, F and G, respectively), and 
have been clearly detected in the Hi map as well. The extent 
of these structures goes well beyond t he galaxy's optical ra- 



dius ( about 20.5 kpc, as estimated by de Vaucouleurs et al 
1991b. While the identification of "E" and "F" is quite clear 



at all SPIRE bands, the detection of structure "G" is more 
problematic, as it is very faint. Nonetheless, the analysis 
of the Hi maps clearly shows an emission which overlaps 
with all three structures, at velocities that are consistent 
with those of Andromeda in that region, making us confi- 
dent also about the presence of the outermost, faintest one. 



Exploiting the velocity information carried by Hi data is 
in this region much simpler, as the south-west side of M31 
has velocities in the opposite direction compared to those 
of the Galactic emission, so that contamination by fore- 
ground Galactic Hi is not an issue. Contours defining these 
three structures are shown in Fig. [8j overlaid to the high 
resolution Hi map. 



In order to quantify in a more objective way the de- 
tection of t hese three struc t ures w e followed the approach 
adopted by Gordon et al. (2009) who, to this aim, used 
the correlation expected between hydrogen column density 
and foreground subtracted FIR maps. We have looked for 
such correlations in regions E, F and G, between the high 
resolution hydrogen map, and SPIRE images. We selected 
these three regions from the higher signal-to-noise Hi map 
based on visually inspecting the morphology of regions with 



Hi column density values greater than ~ 4 x 10 cm" 
(see Fig. |8|. Hence, the exact outline of these regions is, 
to some degree, arbitrary. We present the results for the 
three different regions in Fig. [7j where only SPIRE 350 /im 
is shown as a representative example. Similar results are 
obtained when considering the other two bands, confirming 
our visual identification. The correlation at PACS 160 /im 
is poorer and shows a much larger dispersion and is hence 
not used; at 100 /im we found none. 



While the whole structure of the two rings at 10 and 
15 kpc is visible in all of our maps and can be easily fol- 
lowed all along the five HerscheTs maps, these three newly 
identified structures do not appear to extend all across 
the galaxy, possibly because of their very nature, but also 
due to their lower surface brightness and the high con- 
tamination from Galactic cirrus on the north-east side of 
Andromeda. Despite this, we note that the outermost ring- 
like structure on the northeast side, which is highlighted by 
the 250 /im contour in Fig. [6] is located at roughly the same 
galacto-centric distance as structure E, on the opposite side 
of the galaxy. 



We have also attempted to look for the possible presence 
of dust corresponding with the Hi high-velocity cloud ana- 
lo gues lik ely associated with M31, as identified by Thil ker e"t| 
[ah] ( |2004[ ) . However, after a visual inspection of the SPIRE 
maps, we were unable to identify any obvious signature of 
dust emission corresponding to any of the clouds (see Fig. 
[5]) . As a further check, we have looked at the higher resolu- 
tion WSRT d ata for the M31 H i gh Vel ocity Clouds (HVCs) 
presented in Westmeier et al.] ( |2005[ ), to search for local 
correlations between IR and Hi peaks. Again, we found 
no obvious IR counterparts. On the one hand, we may 
not find any dust emission in these regions possibly be- 
cause of some residual contamination from the foreground 
dust in some regions (especially towards the northeast). 
On the other hand, the foreground can be very variable, 
and the expected emission from these Hi clouds is weaker 
than the foreground. We have indeed checked that the peak 
Hi column density of Davie s' cloud (see Fig. [5]), which is 
the brightest cloud found by Thilker et al.| ([2004 ) , is much 
weaker than the foreground at that position (0.4 x 10 20 



atoms cm and 8.1 x 10 atoms cm , respectively). 
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Fig. 6. Zoom of the north-east edge of M31. The top left panel shows part of the Herschel 250 /im image as a colour 
map, while the same image is shown as contour in the other panels (contours are 3 mJy beam -1 , i.e. ~ 3.8<r, 30 mJy 
beam -1 , and 300 mJy beam -1 ), superimposed to the Hi map at various velocities. The channel velocities are indicated 
in the bottom left corner of each panel. Structures in the IR identified as A to D (see also Fig. [T]and[5| are indicated in 
the panels. The other panels show the SPIRE 250 /jm image as contours and channel maps from the high-resolution Hi 
datacube as colour map. 
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Fig. 7. Pixel-by-pixel correlation between Hi column density, and the flux at 350 /^m for regions E, F and G (see Figs. [I] 
and SI, from left to right, respectively, derived from the foreground-corrected maps. The black lines represent a linear 
fit to the data. 
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Fig. 8. A zoom on the outermost part of the south-west 
region in the high resolution Hi map. Overlapped are the 
contours defining the three outermost regions (see text for 
details on the selection on these regions). 



5. Dust characteristics across the galaxy 

5.1. Dust masses and temperatures 



al. 


2012 


Ford et al. 


in prep. 



performing a more detailed 
analysis and discussion of the dust properties, here we give a 
broad overview of the characteristics of dust in Andromeda, 
both as a function of the galacto-centric distance and con- 
sidering the whole galaxy as well. This approach is similar 
to what has already been done by other studies performed 



Schwering 

TT TT-> <-' 


1987 


Haas et al. 


1998 



2009) 
(see left-hand 

panel of Fig. 9n, of which we report the sizes and character- 
istics in Table jlj the first is a circular one, and it includes 
the innermost dust ring (i? ~ 1.5 kpc). The sizes of the 
other four were obtained by projecting circular apertures 
that define the inner and outer boundaries of the 10 and 



15 kpc rings, plus an outermost aperture including most of 
the structures visible in the SPIRE 250 /im map. Their po- 
sitions were chosen so that their centres coincide with the 
phys ical centres of the s tructures they are meant to identify 
(see Block et al. 2006| ). As a value for the position angle, 
we took that of M31 (38°). As Herschel maps differ both in 
resolution and pixel size, the PSFs in all maps should have 
in principle been matched to the broadest PSF in any map 
(i.e. the 500 ^imPSF). We found this to be not necessary 
as the apertures we have used are much larger than any of 
the pixel and beam sizes. 

We then performed a simple SED fitting assuming that 
the dust thermal emission is well represented by a modified 
black body at a single temperature, with the emissivity 



being described by the following: 
that k = 0.192 m 2 kg" 1 at 350 ^m (|Draine[ [2003 1 



with fcn such 



We could have also used MIPS flux at 70 fim from 
Gordon et al. ( 2006 ) but, as found by other works, this band 
is likely to be more sensitive to warmer dust, and would 
have required the addition of a further modified black body 
component, at higher temperature. Hence, the 70 /im flux 
was only used as an upper limit. We let the temperature be 
a free parameter in our fitting procedure, while the value 
of j3 was kept fixed, and initially assumed to be 2, which 
i s consistent wi t h the value commonly use d in dust models 
praine fc Lee] |1984| |Li fc Drainej |2001[>. The dust mass 
is computed by normalizing the model SED to HerscheFs 
observed datapoints. Uncertainties on the fluxes were con- 
servatively assumed to be 10% for PACS and 7% for SPIRE 
(see ^2.2 and |2.3[ ). In the data fitting, the results of which 
we report in Table [l] we have found the best-fitting temper- 
ature by exploiting a gradient search method. Uncertainties 
have been computed by bootstrapping over 200 SED fits, as 
follows: once the best fit solution is found a set of new, fic- 
titious, datapoints is created by randomly choosing a value 
within the observed errorbars. A best fit model is found 
for these new flux values, and the procedure is repeated 
200 times. The 16% of the best fit models with the lower 
and higher values are discarded, and uncertainties are then 
taken as the differences between the best fit solution and 
the extreme values. 

We compare the value of the dust mass we have ob- 
tained to the one derived by Haas et al. (1998), who use 
a similar approach as ours in their analysis of ISO data. If 
we consider the dust included within the fourth annulus, 
which is the one that better matches their global aperture, 
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Fig. 9. In the left-hand panel the 250 /im image is shown, together with the 5 apertures where fluxes were measured. 
The total aperture corresponds to the outermost elliptical aperture. The best fit values and uncertainties are given in 
Table [TJ In the right-hand panel we show the fit to the observed SEDs as measured in the four apertures: observed 
fluxes (triangles) are superimposed to the modified black body curve (lines) , corresponding to the best fit values of the 
temperature. Different colours and lines correspond to different apertures/annuli. 
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Table 1. Herschel fluxes, as measured within the five annuli we have considered (see Fig. [9J left-hand panel). Dust 
temperature (Tp) and masses (M ), and their uncertainties, are derived from SED fitting (see text for details). We also 
report the mass of hydrogen in the atomic phase, as derived from Hi maps of |Braun et al.| ( |2009[ ) (both corrected an d 
non-corrected for self-opacity), and the mass of gas in the molecular phase derived by CO data (Nieten et al. 2006 1. 
Last two columns give the total gas-to-dust ratios within the given annulus. Results for the largest annuius and for the 
whole galaxy are limited by the size of the CO map (see text for more details): the value indicated with a assumes a 
mass of 5% of total, within that annulus, and the value of the gas-to-dust in the outermost annulus ( b ) was computed 
assuming a negligible contribution from the molecular fraction, and it is hence to be interpreted as a lower limit. Values 
in the columns labeled with * are calculated using the self-opacity corrected maps, "a" and "b" are the major and minor 
axis length, respectively, of the ellipses, given in kpc (first aperture is circular). 



we obtain a total mass of ~ 5.3 x 10 7 M Q which is almost 
a factor of two higher than their 3 x 10 7 Mp) value, and 
also higher than the value derived by |Gordon et al .| (|2006|) 



(4 x 10 7 M Q , at a temperature of 17 K). The temperature 
they obtain from the global SED -16 K- agrees fairly well 
with the value of 16.8 K that we would derive fitting the 
total flux within the fourth annulus, and is somewhat lower 
than t he 18.7 K value found by |Tabatabaei fc Ber khuijscn 
(2010 1 . The difference in the dust mass value might be, at 



least partially, due to the use of different dust e missivity 
coefficients (Haas et al. use the Hildebrand 1983 prescrip- 



tion for calculating the dust mass). Furthermore, the size 
scales used in determining the total dust mass also affect 
the uncertai nty in the total dust mass as shown by |Galliano| 



et al.| ( |2011[ ) to be up to ~ 50%, for example, in the Large 
Magellanic Cloud. The value of this uncertainty will de- 
pend on the morphology of the galaxy and the mixing of 
the different dust phases in the area used for total dust 
mass determination. 

We have also explored different values of /3, namely 2.5, 
1.7 and 1.5. (3 = 2 turned out to be the best suited value for 
reproducing the annuli SEDs. Only if we consider the total 
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aperture, fits with an acceptable \ 2 were obtained also for 
ft = 1.7, yielding in this case a temperature of 17.7 K and 
a total dust mass of 4.7 x 10 7 M . 

The variation of the average temperature -as a function 
of the radius- derived with this analysis is similar, within 
the er rorbars, to that found by Tabatabaei & Bcrkhuijsen 
(2010), derived from MIPS 70 and 160 /xm images, with the 



exception of the central regions in which the temperature 
we found is lower, most likely due both to our coarser sam- 
pling, and to the fact that we are not fitting the MIPS 70 
/im point, which would be more sensitive to warmer dust. 

A s already found in p revious studies (e.g. Haas et al. 

1998[ IGordon et al.| [20061) , the dust emission in the 1.5 kpc 



ring region is dominated by a higher temperature compo- 
nen t. Thi s is also consistent with the study by |Devereux et| 
aD ( 1994 ) , who find that the gas in the innermost regions 



is dominated by an ionized component (seen in Ha emis- 
sion). On the other hand, if we only consider the emission 
from the 10 kpc ring (sampled by our third annulus), we 
obtain an average temperature of 17.0^3 K, and a dust 
mass of 1.95 x 10 7 M Q . Our global SED fitting shows no 
need for an extra, very cold, dust component that, revealing 
its presence only at SPIRE wavelengths, would have been 
impossible to detect by former studies exploiting only data 
up to ~ 170 /im. This issue is tackled in a better detail 



2012), performing a 



in a companion paper (Smith et al. 
pixel-by-pixel analysis. 

5.2. The gas-to-dust ratio 

We have calculated the mass of neutral hydrogen from the 
Hi maps, within the same annuli, and present it in Table [T] 
To take the atomic Helium fraction into account, this value 
is multiplied by a factor of 1.36. 

The CO map was used to derive the mass of molec- 
ular hydrogen by assuming a H2-CO conversion f actor X 
1.9 x 10 20 cm" 2 (K km s" 1 )" 1 (|Strong 



= N(H 2 )/I 


CO = 


& Mattox 


1996 



The CO-to-molecular hydrogen conver- 



rangc 



certainty. Values derived for M31, in particular, can 
from 0.97 x 10 20 to 4.6 x 10 20 cm" 2 (K km s" 1 )" 1 (|Leroy 



et all 120111). In iSmith et all d2012b we have performed a 



detailed, spatially resolved analysis of dust emission con 
sidering all Herschel bands, and we have found values, for 
the Xco factor in the range 1.5-^2.3 x 10 20 , with a best-fit 
value of 1.9 x 10 20 cm- 2 (K km s -1 ) -1 . 

The total molecular gas mass {M mo \ = Mh 2 + Mu e ) 
is then calculated by multiplying the H2 mass value by 
1.36. The CO observati ons only sample th e region which 
includes the 10 kpc ring. |Nieten et al" (2006) estimate that, 
within a region covering the innermost ~ 18 kpc, approx- 
imativcly corresponding to our fourth annulus, a further 
5% of CO mass is expected. In calculating the gas-to-dust 
ratio in the two outermost apertures (and of course in the 
total one), we have applied this correction, also assuming 
that no further molecular gas is expected at larger galacto- 
centric distances. Hence, the values of the gas-to-dust ratio 
for the aperture 4, 5 and the total, that we give in Table 
[l] and that already include this correction, should be con- 
sidered as lower limits, although not too distant from the 
actual values. In fact, it should be noted that the gas mass 
is dominated by the atomic phase, the molecular only con- 
tributing to the total by less than ~ 10%. 



That the gas-to-dust ratio in M31 is not constant 
throughout the galaxy, and increases as a function of the 
radial coordinate, was already shown by several previous 
studies (see e.g. |Bajaja fc Gergelyl |1977[ Walterbos & 



Kennicutt 1988 



JNedialkov et al.l | 2U00p . We estimated an 



average value of ~ 110, calculated over the largest aperture 
(g2d column of Table [l]), and a variation by a factor of ~ 2 
from the innermost ring to the outermost annulus. 

We found a very good correlation between the gas-to- 
dust ratio and the (average) galacto-centric distance, with 
the first quantity exponentially increasing from the centre 
to the outskirts. Even though our analysis is performed 
over very broad apertures, we found a gradient of 0.0163 
dex/kpc, which is con sistent with the value s of the metallic- 



ity gradient found by |Trundle et al.| (| 2002[ ) in M31, ranging 
from -0.027 to -0.013 dex/kpc, under the assumption that 
the amount of metals locked up into dust grains is a con- 



stant fraction of the mass of interstellar metals ( Edmunds 
|2001[ ). Should this be the case, then the metallicity gradient 
would be the opposite as the one for the gas-to-dust. 



Braun et al. (20091 also provide maps where hydrogen 
self-opacity is taken into account. It is quite clear that the 
simple 21 cm integral is a lower limit to the Hi mass, and 
M31 is one of the few cases where we have the data to ac- 



count for this problem (Braun 2012). In fact, had we taken 



this into account, the atomic mass would have increased by 
a factor ranging from ~ 5% in the innermost regions, to 
~ 27% in the most distant annulus. In Table [l] we report 
the Hi mass and the gas-to-dust ratio calculated from the 
self-opacity corrected maps (columns indicated with a *). 



5.3. Dust in the outskirts 

We will now briefly characterize the properties of dust that 
is detected in the outermost regions. In Sec. |4.1| and 472] we 
have identified structures in the outskirts of M31 disk, that 
are likely to belong to Andromeda. As even after subtract- 
ing a model for the Galactic cirrus emission the situation 
remains confused in the North-East side, we will focus on 
the structures that are identified on the southern part, as 
they are much less affected by foreground dust emission and 
which we have labelled as E, F and G. 

In Sec. |4.2| we have shown that there is a significant 
correlation between the hydrogen column density and FIR 
emission at SPIRE wavelengths for these arc-like shaped 
structures. Using the same technique that is described in 
Sec. [5j we fit the observed fluxes measured within the ar- 
eas defined by these regions (see Fig. [8]), by means of a 
modified black body SED. As there is no clear detection 
at PACS wavelengths, we are forced to limit ourselves to 
SPIRE wavelengths. 

As before, we have considered four values for the emis- 
sivity spectral index ft, namely 1.5, 1.7, 2.0 and 2.5. 
Acceptable fits were achieved with all the four values of 
ft for region E, with ft = 1.5 being slightly preferred. The 
SEDs of regions F and G turned out to be compatible with 
ft = 1.5 and 1.7 only, with the first value yielding, again, 
slightly better fits. Considering only ft — 1.5, we found dust 
masses in these three regions of 1.23 x 10 6 M , 5.89 x 10 5 
M Q and 3.55 x 10 5 M , respectively, and temperatures of 
14.4, 13.1 and 11.6 K. Measuring the hydrogen mass on the 
same regions, we derive Hi-to-dust ratios of 202, 204 and 
274, respectively. 
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6. Summary and conclusions 

We have presented and described Herschel PACS (100 and 
160 /zni) and SPIRE (250, 350 and 500 /zm) observations of 
a ~ 5.5° x 2.5° field centered on the local galaxy M31. Using 
high spatial and velocity resolution Hi maps, roughly cov- 
ering the same area, we have identified the dusty structures 



~ 3° x 1° reg 


on (|Haas et al.| 1998 Gordon et al. 


200G 


Barmby et ah] 


200G: 


Krause et al. in prep.p. 



We have performed a broad analysis of the characteris 
tics of the dust as a function of the radius. We have consid- 
ered 5 main regions, described by projected circular aper- 
tures (plus an elliptical one), sampling the most significant 
morphological structures: the inner 1.5 kpc ring, the intra- 
arm region enclosed within the 10-kpc ring, the 10-kpc 
ring itself, the outer (~ 15 kpc) ring, and a fifth one, con- 
taining the outermost regions. A simple one-temperature 
modified black body model is able to successfully repro- 
duce the observed PACS/SPIRE SED. Average quantities, 
such as the dust temperature, and the gas-to-dust ratio, 
are derived, and they are found to strongly depend on the 
radius. The gas-to-dust ratio, in particular, is changing by 
a factor of ~ 2, following an exponential increase, from the 
inner regions to the outskirts, as expected from the pres- 
ence of a metallicity gradient with the opposite sign, across 
the galaxy. 

According to our simple emission model the two main 
rings, at 10 and 15 kpc from the centre, contain the bulk 
(about 78%) of the dust in all the galaxy. The dust sur- 
face density is the highest in within the 10 kpc one, 0.41 
M Q /pc 2 , being higher by a factor of ~ 4 with respect to 
the average value, and also higher than in the 15 kpc ring, 
where it is about 0.25 M /pc 2 . 

The dust temperature is highest in the innermost aper- 
ture, but the value we derive from our SED fitt ing, 22 K, is 



considerably lower than previous stud ies (e.g. Haas et al. 
1998; Tabata baei fc Ber khuijsen| |2010[ derive temperatures 
of ~ 28 and ~ 30 K, respectively, for the nuclear dust). This 
is most likely due to the fact that, ignoring data shortwards 
100 /im, we are not sensitive to the presence of a warmer 
component (e.g. Bendo et al. 2010): this would, anyway, 



constitute a minor fraction of the total dust content. 

We will give a more detailed description of the I R mor- 



1 



m prep. 



, where 



Herschel maps 



phology in a forthcoming paper (Kirk et al. 
an analysis of the structures on deprojectec 
will be performed. 

The maps we describe in this work present two distinct 
situations: in the north-eastern part, the combination of 
Andromeda's systemic and rotation velocities, yielding val- 
ues that are similar to those expected for the interstellar 
medium in our galaxy, together with the presence of a sub- 
stantial amount of gas and dust, makes the identification of 
emission regions that actually belong to Andromeda very 
difficult, and sometimes ambiguous. On the other hand, at 
the southern side, not only is the velocity of gas and dust 
that belong to Andromeda expected to have a clearly differ- 
ent pattern with respect to that of our Galaxy, but is also 
much less affected by the presence of interstellar medium 
(see e.g. Fig. [5]for a panoramic view of the M31 neighbour- 
hood). 

In this region, we identify in our maps three different 
ring-shaped, concentric structures - which we label E, F 
and G - that can be safely attributed to M31, based on the 



shape and velocity of their Hi counterparts, extending out 
to a radius of ~ 31 kpc. Dust in the outskirt of Andromeda 
had been already indirectly detected, through extinction 
of b oth s tars and background galaxies, by Cuilland re et| 



al 



(2001), who exploited deep optical photometry in a 
28' x 28' field roughly centred on our "E" region. This 
res ult, re cently confirmed on a smaller scale by |Bernard et 
al. ( 2012[ ) (their "Outer Disc field"), implies an atomic gas- 
to-dust ratio about three times higher with respect to the 
average Galactic value, and is considered as an upper limit 
since it only takes into account the dust in front of the ob- 
served stars. Instead, the Hi-to-dust ratio we estimate for 
this particular region, ~ 204, is a factor of ~ 1.5 larg er than 
the 137 valu e estimated for the MW (see Table 2 in Draine 

y solar 



et al. 



met a 



2007 1 , that one would expect given the roughl; 
icity of stars observed in a nearby field. This value 
can provide a more stringent constraint to the amount of 
dust in this region, and to the chemical enrichment of the 
interstellar medium in Andromeda's outskirts. 
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